Charge Order and the Origin of Giant Magnetocapacitance in LuFe 2 04 
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The nature of the charge order in the charge frustrated compound LuFe2C>4 and its effect on 
magnetocapacitance were examined on the basis of first-principles electronic structure calculations 
and Monte Carlo simulations of electrostatic energy. Our work shows that two different types of 
charge order of almost equal stability (i.e., y3 x y/i and chain types) occur in the Fe2C>4 layers of 
LuFe204, and that the ground state of LuFe204 has a ferrielectric arrangement of the Fe204 layers 
with y/3 x \/3 charge order. The giant magnetocapacitance effect of LuFe204 at room temperature 
is accounted for in terms of charge fluctuations arising from the interconversion between the two 
types of charge order, that becomes hindered by an applied magnetic field. 

PACS numbers: 71.20.-b,71.45.Lr,77.80.-e,77.84.-s 
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Ferroelectric (FE) oxides are essential components in 
a large number of applications [H . In traditional FE ma- 
terials like BaTi03, the ferroelectricity is driven by the 
hybridization of the empty d orbitals of Ti 4+ with the 
occupied p orbitals of the oxygen anions @. Recently, 
a mixed- valence compound LuFe204 with the average 
valence Fe 2,5+ was found to exhibit ferroelectricity as- 
sociated with the charge order (CO) leading to Fe 2+ and 
Fe 3+ ions Q . Subramanian et al. reported that at room 
temperature (RT) the dielectric constant of LuFe2C"4 de- 
creases sharply when a small magnetic field is applied [1] . 
This suggests a strong coupling between spin moment 
and electric dipole at RT, and hence potential applica- 
tions of LuFe204 in which the charge and spin degrees of 
freedom of electrons can be controlled. 

At RT LuFe204 has a hexagonal layered structure 
(space group R3m, a = 3.44 A, and c = 25.28 A) in 
which all Fe sites are crystallographically equivalent 5]. 
LuFe204 is an insulator [6j], and undergoes a three- 
dimensional CO below 330 K [1, 0, S| as well as a two- 
dimensional ferrimagnetic order below 240 K |. In 
LuFe2 04, layers of composition Fe204 alternate with lay- 
ers of Lu 3+ ions, and there are three Fe204 layers per unit 
cell (Fig. 1(a)). Each Fe2 04 layer (referred to as the 
W-layer) is made up of two triangular sheets of corner- 
sharing FeOs trigonal bipyramids (Fig. l(b)-(d)). 

LuFe204 exhibits apparently puzzling physical prop- 
erties, and the nature of its CO is not unequivocal. It 
is generally believed that lattice distortions accompany 
a CO [9(. Thus, it is unclear why LuFe204 is insulat- 
ing above RT despite that LuFe204 adopts the structure 
in which all Fe sites are equivalent in this temperature 
region. The giant magnetocapacitance effect found for 
LuFe204 can be understood by supposing that the charge 
fluctuation (CF) of LuFe204 arising from its charge frus- 
tration is sharply reduced by an external magnetic field. 
However, it is unclear by what mechanism this happens. 
As to the nature of the CO in LuFe204, there is a con- 
troversy. By analogy with the stable long-range order of 
Ising spins in a triangular lattice antiferromagnet (TLA) , 



Yamada et al. proposed a model CO structure in which 
the two triangular sheets of a W-layer do not have the 
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same number of Fe 2+ and Fe 3+ ions, i.e., [Fe 2+ ]:[Fe 
1:2 in one triangular sheet, and 2:1 in another triangular 
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FIG. 1: (Color online) (a) Hexagonal unit cell structure of 
LuFe2 04, where grey, white and red circles represent the Lu, 
Fe and O atoms, respectively. Schematic representations of 
the (b) CO-I, (c) CO-II and (d) CO-III structures, where the 
small and large circles refer to the Fe ions in the upper and 
lower triangular sheets of a W-layer, respectively. The Fe 2+ 
and Fe 3+ ions are represented by empty and filled circles, 
respectively. Schematic representations of the (e) FE and (f) 
FIE arrangements of the dipoles of the W-layers in the CO-I 
structure. 
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sheet (i.e., the CO-I structure in Fig. 1(b)) [U$ . The 
y/3 x V3 superstructure of this CO model is compatible 
with the experimental observation. Subramanian et al. 
proposed a different CO structure in which one triangu- 
lar sheet of a W-layer has only Fe 2+ ions and the other 
triangular sheet has only Fe 3+ ions (i.e., the CO-II struc- 
ture in Fig. 1(c)) 4], on the basis of the fact that the 
closest Fe-Fe distances in LuFe2 04 occur between adja- 
cent triangular sheets of a W-layer rather than within 
each triangular sheet. 

In this Letter, we probe the nature of the CO and 
the origin of the magnetocapacitance in LuFe204 on the 
basis of first principles electronic structure calculations 
and Monte Carlo (MC) simulations of electrostatic in- 
teractions. Our study shows that the ground state of 
LuFe204 has the CO-I structure, the CO-II structure is 
unstable, and another CO structure different from CO-I 
and CO-II is very close in energy to the CO-I structure. 
The presence of two different CO structures close in en- 
ergy is found crucial for the magnetocapacitance effect of 
LuFe 2 4 . 

Our spin-polarized density functional theory calcu- 
lations were performed on the basis of the frozen- 
core projector augmented wave method [111 encoded in 
the Vienna ab initio simulation package [12| using the 
generalized-gradient approximation (GGA) [131 ] and the 
plane-wave cut-off energy of 400 eV. To properly de- 
scribe the strong electron correlation in the 3d transition- 
metal oxide, the GGA plus on-site repulsion U method 
(GGA+U) was employed 14| with the effective U value 
(U e ff = U — J ) of 4.61 eV. Calculations with various 
U e f f values show that our main results remain valid when 
U e ff is varied between ~3.6 and ~5.6 eV. In the follow- 
ing we report only those results based on C/ e // = 4.61 eV. 
To simplify our discussion and reduce the computational 
task, we explore the CO structures of LuFe2 04 under 
the assumption that the spins of LuFe204 have a ferro- 
magnetic (FM) ordering. This assumption is reasonable 
because the energy scale associated with different spin 
arrangements is much smaller than that associated with 
different CO's. An additional restriction of our calcu- 
lations is the use of the experimental lattice constants. 
Our full geometry optimization of LuFe 2 04 with GGA 
calculations leads to the lattice constants that are very 
close to the experimental values. 

The GGA+U calculation for the non-CO state of 
LuFe2 04 with the room-temperature crystal structure 
and hence the hexagonal symmetry shows that all Fe ions 
are in the high spin state. The plots of density of states 
(DOS) calculated for the non-CO state (Fig. 2(a)) pre- 
dict a metallic behavior for LuFe2 04, in disagreement 
with experiment. This is not surprising because the Fe 
ions are in the valence state of 2.5+ in the absence of CO. 
We found it impossible to produce an insulating band gap 
for any CO structure as long as the structure keeps the 
hexagonal symmetry. For a transition metal ion in a trig- 



onal bipyramidal crystal field, the five d states are split 
in three groups, i.e., {d xy , d x 2_ y 2}, {d xz , d yz } and {d z s}. 
The partial DOS plots shown in Fig. 2(b) indicates that 
the down-spin bands leading to a metallic character arise 
from the d xy and d x 



orbitals, i.e., the {d xy , d x 2_ y 2} 



states have a lower energy than do the {d xz , d yz } states 
for the FeOs trigonal bipyramids in LuFe2 04. 
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FIG. 2: (Color online) DOS of LuFe204 calculated for the 
RT structure without CO. The DOS was calculated with 0.1 
eV broadening. 

To search for possible low-energy CO patterns, we re- 
sort to the classical MC simulation method by consider- 
ing only the intcrsitc Coulomb repulsion. This simplified 
approach is reasonable for the purpose of finding stable 
CO patterns. As the energy reference of this calcula- 
tion, each Fe 25+ site of the non-CO state is assumed to 
carry zero charge, so that, after a CO takes place, Fe 3+ 
and Fe 2+ sites will carry 0.5<5 (0 < 5 < 1) and —0.5(5 
charges, respectively (6 denotes the degree of the charge 
transfer). For simplicity, we assume that the charge on 
oxygen is independent of i5. We perform MC simulations 
using a 2^/3 x 2\/3 supercell with periodic boundary con- 
ditions and evaluate the long-range Coulomb interaction 
using the Ewald sum method. There occur two differ- 
ent low-energy CO patterns. The most stable one is 
a chain-like CO structure (hereafter the CO-III struc- 
ture) with energy — 1.6015 2 eV per formula unit (FU), in 
which chains of Fe 2+ ions alternate with those of Fe 3+ 
ions in each triangular sheet (Fig. 1(d)). Given the trig- 
onal symmetry of each W-layer, there are three differ- 
ent, but equivalent, ways of choosing the chain direction 
in each W-layer. When the three different chain orien- 
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tations occur randomly, the overall superlattice diffrac- 
tion pattern will exhibit a y/3 x V3 structure. The COT 
structure proposed by Yamada et al. is found to have a 
slightly higher energy, i.e., -1.3796 2 eV/FU. (The sta- 
bility of this CO cannot be deduced by analogy with the 
V3 x \/3 spin configuration of Ising spins that occurs in 
a TLA when the nearest-neighbor (NN) spin exchange is 
antiferromagnetic and the next-nearest-neighbor (NNN) 
spin exchange is FM. This spin configuration becomes the 
ground state only when the NNN interactions are stabi- 
lizing. In the CO-I structure, the NNN interactions are 
destabilizing due to Coulomb repulsion.) Our MC simu- 
lations show that the CO-II structure proposed by Sub- 
ramanian et al. is found to be highly unstable, namely, 
its energy is even higher than the non-CO structure by 
0.496.5 2 eV/FU. The instability of the CO-II structure is 
due to the large positive electrostatic energy between the 
ions of the same charge in each triangular sheet. We also 
consider MC simulations using another model that takes 
into consideration only the NN electrostatic interactions 
within each triangular sheet and between adjacent trian- 
gular sheets of every W-layer. Qualitatively, this model 
leads the same results as described above. 

To more accurately probe the stabilities of the CO-I, 
CO-II and CO-III structures with respect to the non- 
CO structure, we carry out GGA+U calculations for 
LuFe2 04 on the basis of its crystal structure determined 
at RT. Our GGA+U calculations for the CO-II structure, 
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FIG. 3: (Color online) DOS of LuFe 2 4 calculated for (a) the 
relaxed CO-I and (b) the relaxed CO-III structures. The blue 
regions highlight the PDOS of the Fe 2+ ion. The inset shows 
the charge density plot calculated for the occupied down-spin 
d-states. The DOS was calculated with 0.1 eV broadening. 



with an initial guess of the charge density expected for 
it, always converged to the non-CO structure. Therefore, 
it is concluded that the CO-II structure is not stable, as 
found from our MC simulation. Since there are three W- 
layers per unit cell and since each W-layer has a nonzero 
dipole moment in the CO-I structure, there are two differ- 
ent arrangements of the W-layers in the CO-I structure, 
i.e., the FE arrangement (Fig. 1(e)) and the ferrielectric 
(FIE) arrangement (Fig. 1(f)). For the FE CO-I and the 
CO-III structures, our GGA+U calculations were carried 
out in two steps. In the first step, the FE CO-I or the 
CO-III structure was introduced in LuFe204 without al- 
lowing the crystal structure to relax. Such a CO induced 
solely by electrostatic interactions was first proposed by 
Attfield et al. 11511 and later computationally realized by 
Leonov et al. [16| in their studies of the CO phenomenon 
in Fe20B03. In the second step, the crystal structure 
with the FE CO-I or the CO-III structure was completely 
optimized. Thus, the energy gain AE of the FE CO-I or 
the CO-III structure relative to the non-CO structure is 
written as AE = AE± + AE 2 , where AEi and AE 2 are 
the energy gains obtained in the first and second steps, 
respectively. A-Ei arises from electrostatic interactions, 
and AE 2 from the geometry relaxation. 

Our GGA+U calculations of the first step show that 
the FE CO-I and the CO-III structures are stable in the 
absence of geometry relaxation [3] (with AE\ = 239 and 
219 meV/FU, respectively), and are an insulator (with 
band gap of 0.18 and 0.15 eV, respectively). The forces 
acting on the atoms calculated for both CO structures 
with no geometry relaxation are large hence indicating 
instability of the "frozen" structure. Our GGA+U cal- 
culations of the second step show that the AE2 values 
of the FE CO-I and the CO-III structures are 163 and 
165 meV/FU, respectively. Thus, after structural relax- 
ation, the FE CO-I structure remains only slightly more 
stable than the CO-III structure (by 19 meV/FU). An- 
other important finding of our calculations is that the 
energy gain resulting solely from electrostatic interac- 
tions is greater than that from geometry relaxation (i.e., 
AEi > AE 2 ), in contrast to the case of Fe 2 OB0 3 , for 
which our calculations showed that AEi <C AE 2 (i.e., 77 
vs. 272 meV) !17j- Bond- valence-sum calculations [3] 
for the optimized structures of the FE CO-I and the CO- 
III structures show that the valence states of the "Fe 2+ " 
and "Fe 3+ " sites are close to the nominal +2 and +3, 
respectively. Both CO structures have a larger band gap 
after geometry optimization (i.e., 0.70 and 0.54 eV for 
the FE CO-I and the CO-III structures, respectively). 
The DOS plots calculated for the FE CO-I and the CO- 
III structures with the relaxed structures are presented 
in Fig. 3. All the occupied up-spin 3d bands of the 
Fe 2+ ions lead to a spherical charge distribution, and so 
do those of the Fe 3+ ions. It is the occupied down-spin 
3d bands of the Fe 2+ ions that provide an anisotropic 
charge distribution arising from the Fe 3d orbitals and 
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hence information about the orbital order (00). This is 
shown in the insets of Fig. 3, which reveal that the two 
triangular sheets of a W-layer have an identical 00 in 
the CO-III structure, but different OO's in the FE CO-I 
structure. In the FE CO-I structure, the orbitals of the 
{d xy , d x 2_ y 2} type are involved in the 00 of the trian- 
gular sheet of containing fewer Fe 2+ ions than Fe 3+ ions, 
but those of the {d xy , d x 2_ y 2} and {d xz , d yz } types in 
the other triangular sheet. In the CO-III structure, the 
orbitals of the {d xy , d x 2_ y 2~\ and {d xz , d yz } types are 
involved in the 00 of both triangular sheets of a W- 
layer. The OO's found for the unrelaxed structure are 
quite similar to those found for the optimized structures 
described above. 

Experimentally, a large spontaneous electric polariza- 
tion (EP) was found in LuFe204 [3J. Of the two stable 
CO structures, CO-I and CO-III, only the CO-I structure 
has a nonzero EP. For the calculations of the spontaneous 
EP of the FE CO-I structure, we use the Berry phase 
method [23]. The calculated EP should be compared 
with the experimental one measured at a low tempera- 
ture where there is no CF. The spontaneous EP along 
the c-axis is calculated to be 52.7 /iC/cm 2 , which is 
much greater than the experimental value of about 25 
liC/cm 2 measured at 77 K [3|. To resolve this discrep- 
ancy, we consider the FIE CO-I structure (Fig. 1(f)). 
Our GGA+U calculations with complete geometry opti- 
mization show that the FIE CO-I structure is more stable 
than the FE CO-I structure by 12 meV/FU. The spon- 
taneous EP along the c-axis calculated for the FIE CO-I 
structure is 26.3 /iC '/ 'cm 2 , in excellent agreement with ex- 
periment. Thus, from the viewpoint of the total energy 
and the spontaneous EP, it is concluded that LuFe2 04 
has the FIE CO-I structure. Note that the spontaneous 
EP for the FIE CO-I structure is smaller than that of FE 
CO-I structure by a factor of approximately two, instead 
of three expected from Fig. 1(e), (f). This is due to the 
fact the two crystal structures, which are each separately 
optimized, are slightly different. 

Now we turn attention to the probable origin of the 
giant magnetocapacitancc effect of LuFe204 at RT. Ac- 
cording to our calculations, the two different CO struc- 
tures, CO-I and CO-III, are very close in energy. Given 
the room-temperature structure of hexagonal symmetry, 
the CO-III state is higher in energy than the FE CO- 
I state only by 20 meV/FU. Thus, at RT, there should 
occur CFs associated with the interconversion between 
the two different CO states, which should form different 
domains separated by domain boundaries. Then, due to 
the large polarizability caused by the CFs, the dielec- 
tric constant of LuFe2 04 should be very large. Under an 
external magnetic field the Zeeman effect should prefer- 
entially stabilize one of the two CO states because the 
two states are most likely to have different total spin 
moments. Consequently, an external magnetic field will 
reduce the extent of CF and hence decrease the electron 



polarizability. The a.c. dielectric dispersion observed for 
LuFe204 can be understood in terms of the dielectric 
response and the motion of the FE domain boundaries 
between the CO-I and CO-III states [3]. At present we 
cannot answer the question of which CO state will be 
preferentially stabilized by an external magnetic field, 
because only the FM spin arrangement was considered 
in the present work. Further studies are necessary to 
address this question. 

We note that the spontaneous EP of LuFe2 04 is much 
greater than that of the recently discovered multiferroic 
materials (such as TbMn03 [___[), so LuFe2 04 represents 
a promising candidate for novel magnetoelectric devices. 
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